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ABSTRACT 

The effect of the Hall term on the evolution of the magnetorotational instability 
(MRI) in weakly ionized accretion disks is investigated using local axisymmetric 
simulations. First, we show that the Hall term has important effects on the MRI when 
the temperature and density in the disk is below a few thousand K and between 10 13 
and 10 18 cm -3 respectively. Such conditions can occur in the quiescent phase of dwarf 
nova disks, or in the inner part (inside 10 - 100 AU) of protoplanetary disks. When 
the Hall term is important, the properties of the MRI are dependent on the direction 
of the magnetic field with respect to the angular velocity vector 17. If the disk is 
threaded by a uniform vertical field oriented in the same sense as 17, the axisymmetric 
evolution of the MRI is an exponentially growing two-channel flow without saturation. 
When the field is oppositely directed to 17, however, small scale fluctuations prevent 
the nonlinear growth of the channel flow and the MRI evolves into MHD turbulence. 
These results are anticipated from the characteristics of the linear dispersion relation. 
In axisymmetry on a field with zero-net flux, the evolution of the MRI is independent 
of the size of the Hall term relative to the inductive term. The evolution in this case is 
determined mostly by the effect of ohmic dissipation. 

Subject headings: accretion, accretion disks — diffusion — instabilities — MHD - 
turbulence 

1. INTRODUCTION 

The structure and evolution of accretion disks are largely determined by angular momentum 
transport processes. One of the most promising processes is MHD turbulence driven by the 
magnetorotational instability (MRI) (Balbus & Hawley 1991). In ideal MHD, the growth rate of 
the MRI is of the order of the orbital frequency O, and the characteristic wavelength of the most 
unstable mode is 2ttva/Q, where v\ is the Alfven speed. The nonlinear regime of the MRI has 
been well studied in ideal MHD using numerical simulations. However, in some systems, accretion 
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disks are expected to be only partially ionized, and non-ideal MHD effects, which generally 
suppress the growth of the MRI, must be considered. For example, the low temperatures of 
accretion (protoplanetary) disks around young stellar objects make thermal ionization processes 
ineffective, so that the abundance of charged particles is very small, and non-ideal MHD effects are 
important (Gammie 1996; Stone et al. 2000). Another example is provided by dwarf nova systems 
in quiescence. The temperature of the disk in this case can be well below 10 4 K and the disk is 
only weakly ionized, so that non-ideal MHD effects are again important (Gammie & Menou 1998). 

There are three regimes in non-ideal MHD associated with the relative importance of different 
terms in the generalized Ohm's law (see §2); they are the ambipolar diffusion, ohmic dissipation, 
and Hall regimes. Which term dominates depends on the ionization fraction and density of the 
gas. Ambipolar diffusion is most important in regions of relatively low density and high ionization 
(e.g., Regos 1997). The linear properties of the MRI in the ambipolar regime have been explored 
by Blaes k, Balbus (1994); they find unstable modes exist when the collision frequency of an ion 
with neutrals is higher than the orbital frequency. The nonlinear evolution of the MRI in this 
regime was examined by Hawley k, Stone (1998) using two fluid simulations. They found that 
when the coupling between ions and neutrals is weak, the turbulence in the ionized component of 
the plasma excited by the MRI does not affect the motion of neutrals very much, thus significant 
angular momentum transport requires a greater coupling between the ions and neutrals than that 
required for linear instability. Brandenburg et al. (1995) and Mac Low et al. (1995) also studied 
the effect of the ambipolar diffusion in the strong coupling limit in a few models. 

Ohmic dissipation becomes important when the ionization fraction of the gas is very low. 
In this case, a linear analysis (Jin 1996; Sano & Miyama 1999) shows that small wavelength 
perturbations are damped, and the characteristic wavelength of the MRI increases in comparison 
to the ideal MHD case. The axisymmetric 2D evolution of the MRI demonstrates that nonlinear 
saturation can occur due to ohmic dissipation (Sano, Inutsuka, & Miyama 1998), even though the 
corresponding ideal MHD cases show an ever-growing channel flow without saturation (Hawley & 
Balbus 1992). Fleming, Stone, & Hawley (2000) examined the nonlinear evolution using local 3D 
simulations; they found that dissipation weakens the MHD turbulence. For significant turbulence 
and angular momentum transport to occur, a critical value for the magnetic Reynolds number 
must be exceeded, and this value depends on the field geometry in the disk. 

Recently, linear analyses of the MRI in the Hall regime have been presented by Wardle (1999) 
and Balbus & Terquem (2001). The maximum growth rate and characteristic wavelength of the 
MRI are strongly modified by the Hall effect. Most interesting is that the linear properties of the 
instability depend on the direction of the magnetic field. This is because the dispersion relation 
for incompressible Alfven waves traveling along field lines is quite different in Hall MHD. In 
particular, the left- and right-circularly polarized Alfven waves have different phase velocities, and 
these two waves interact with the Coriolis force in the disk in different ways. One of these two 
waves is commonly referred to as the whistler wave. 



-3- 



The purpose of this paper is to investigate the effect of the Hall term on the nonlinear 
evolution of the MRI in axisymmetry. When the Hall term is important, ohmic dissipation often 
cannot be neglected (Balbus & Terquem 2001). Thus, we solve an induction equation that includes 
both ohmic dissipation and the Hall effect. The plan of this paper is as follows. We examine when 
the Hall term becomes important in dwarf nova and protoplanetary disks in §2. Our numerical 
method and the initial conditions are described in §3. The results of 2D MHD simulations are 
presented in §4 for both a uniform and a zero-net flux vertical field. The application of the results 
to actual accretion disks is discussed in §5, and our results are summarized in §6. 



2. THE HALL REGIME IN ACCRETION DISKS 

In this section we examine when the Hall term becomes important in actual accretion disks. 
We consider a weakly ionized fluid composed of ions, electrons, and neutrals. Charged dust grains 
can be important in some situations (Wardle & Ng 1999; Sano et al. 2000), because the ionization 
fraction of the plasma is strongly affected by the abundance and size distribution of the grains 
through recombination processes on grain surface. If dust grains with the interstellar MRN size 
distribution (Mathis, Rumpl, & Nordsieck 1977) are considered, the Hall term is important for 
densities between 10 7 and 10 11 cm~ 3 (Wardle & Ng 1999). At these densities, negatively charged 
grains and positively charged ions are the dominant charge carriers. At higher densities than 10 11 
cm -3 , both the negative and positive charge carriers are grains, and the dynamics is dominated 
by ohmic dissipation. The expected density of protoplanetary disks is higher than the Hall regime 
in a dusty plasma, thus if dust is present most of the disk will be in the ohmic dissipation regime 
(Wardle & Ng 1999). However, this assumes the dust grains are well mixed, and are unevolved. 
In fact, dust grains in the disk may grow in size through collisions and/or sediment toward the 
midplane. These evolutionary effects reduce the abundance of charged grains, so that ions and 
electrons become the dominant charge carriers even in very dense regions (Sano et al. 2000). On 
the other hand, whether the gas flow in the disk is laminar or turbulent can largely affect the 
evolution of the grains (e.g., Cuzzi, Dobrovolskis, & Champney 1993; Hodgson & Brandenburg 
1998). Self-consistent models of dusty disks which include the effects of gravitational settling, 
turbulent mixing, and grain growth are beyond the scope of this paper. The analysis presented 
here, which does not include the effect of dust grains, could be applicable to protoplanetary disks 
in the late stages of evolution. Dust is not expected to be important in dwarf nova disks. 

Charge neutrality, n e = rii, is assumed, where n e and rii are the number densities of electrons 
and ions respectively. The coupling between charged particles and neutrals depends on the collision 
rate (av) a , where a is i (ions) and e (electrons). Here, we adopt the values (av)i = 1.9 x 10~ 9 cm 3 
s- 1 and (av) e = 8.28 x lQ-^T 1 / 2 cm 3 s" 1 (Draine, Roberge, & Dalgarno 1983). The imperfect 
coupling between electrons and neutrals results in a finite electrical conductivity, 
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where e and m e are the electron charge and mass. The electron- neutral collision frequency is 
v e — n n(o'v) e , where n n is the number density of neutrals. The associated magnetic diffusivity is 

c 2 
Aira c 

where c is the speed of light. The cyclotron frequencies of the ions and electrons are 

- eB 

m a c 

where a is i and e. 

The induction equation including the terms of the ohmic dissipation, the Hall effect, and the 
ambipolar diffusion is given by 



47T77J J X B (J X B) X B 

v x B + v 



dB „ 

at c en e ciPiP 

where v is the neutral velocity and 

.7 = 

47T 



(4) 



J = .£. (V x B) (5) 



is the current density. The drag coefficient is 

7 = ^=^, (6) 
rrii + m n 

where rrii and m n = fimn are the ion and neutral particle mass, /i is the mean molecular weight of 
the gas, and mn is the mass of a hydrogen atom. 

Here we estimate the relative magnitude of the four terms on the right hand side of the 
induction equation. Following Balbus & Terquem (2001), we denote these terms (reading from left 
to right) as I (Inductive), O (Ohmic), H (Hall), and A (Ambipolar). The ratios A/H and H/O 
are given by 

A- w ci /„n 

77 = — . (7) 
H 7p 

and 

H w ce , . 

O v e 

Note that these ratios are determined only by micro-physical quantities. In this paper we define 
the Hall regime as the conditions under which the Hall term is the largest of the three non-ideal 
MHD effects. From equations (|?|) and (|8|), in the Hall regime electrons are coupled to the magnetic 
field {w ce /v e > 1), but ions are not (w C i/jp < 1) because of the collisions with neutrals. In low 
density regions, both ions and electrons are coupled to the magnetic field: this is the ambipolar 
diffusion regime (A > H > O). On the other hand, in the high density regions, both the electrons 
and ions are decoupled from the field, leading to the ohmic dissipation regime (O > H > A). The 
density range of the Hall regime is therefore intermediate between these two (see below). 
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The importance of each non-ideal MHD effect relative to the inductive term is 

O rj rjO, 

H cB cBn X 



. - (10) 
/ 4iren e VL Aixen e vj^ 2 

and 

A _ B 2 n 
I AnjpipVL 7 pi 

for the ohmic dissipation, the Hall effect, and the ambipolar diffusion, respectively, where V and L 
are typical values for the velocity and length scale in the fluid. We assume that V is of the order 
of va, and that L = va/Q, i.e., the characteristic scale of the MRI. Equations (|9|) and (10) give 
the definition of two important non-dimensional parameters; the magnetic Reynolds number Reu 
and the Hall parameter X. The relation between these two parameters is given by 

H _w ce _ XRe M h0 , 
O ~ v e ~ 2 ' [U) 

This is independent of the ionization fraction and the choice of the typical scales V and L, and 
constrains the values of X and Reu in the Hall regime. Note that the definition of the Hall 
parameter X is the same as x in Balbus & Terquem (2001) (in this paper we shall use x to 
denote one of the coordinates in our computational domain). The importance of the Hall effect 
is estimated by the ratio H/I, or X. In order to evaluate this ratio in actual disks, we need to 
know the number density of electrons n e . In the following subsections, we calculate the ionization 
fraction, n e /n n , in dwarf nova and protoplanetary disks. 



2.1. Dwarf Nova Disks in Queiscence 

In outburst, the typical temperature of dwarf nova disks is over 10 4 K so that the gas is fully 
ionized (e.g., Cannizzo 1993). However, in quiescence the temperature is only a few thousand K. 
The main source of free electrons at this temperature is thermal ionization of Na, Al, and K. 
Assuming the mean mass of ions is mj = 30mn, and the neutrals are hydrogen atoms (p = 1.27), 
equations (|7j) - (^) give 

and 

O nfn e \- l f T \-^ 2 ( c s \ 2 ( M X 1 / 2 / r \ ~ 3 / 2 , x 

7 = " x 10 " (t) Kwk) i£ {Mi) (10^=) < 15 > 
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where r is the distance from the central star, which has mass M. For M = M@, T = 10 3 K, and 
va/cs = 1) the density range of the Hall regime is 10 13 <^ n n <^ 10 18 cnr 3 . The typical density of 
a dwarf nova disk in quiescence is n n ~ 10 18 cm -3 (Gammie & Menou 1998), so that the disk is in 
the Hall regime. The ratio H/O varies from unity to 160(T/10 3 K) -1 / 2 in the Hall regime, so that 

2 £ \X\Re M £ 320 [j^j • (16) 

The ratios of non-ideal MHD terms to the inductive term in the induction equation (Q) are 
functions of the ionization fraction n e /n n . The electron density n e in dwarf nova disks comes from 
solving the Saha equation for Na, Al, and K assuming solar abundances. Figure |l| shows the ratios 
H/I, O/I, and A/I as a function of the number density of neutrals n n . We assume c s /va = 1 and 
r = 10 10 cm in this figure. When the ratio H/I takes a value greater than unity, the Hall effect 
changes the linear character of the MRI (Balbus & Terquem 2001). This ratio is found to be very 
sensitive to the temperature of the disk. 

For T = 3000 K, the typical value of |X| in the Hall regime is 10 -3 at n n ~ 10 18 cm~ 3 , so 
that the Hall term is dominant, but not large enough to affect the linear growth of the MRI. The 
magnetic Reynolds number R&m = (Q/I)~ l ~ 10 4 in this case, implying that ohmic dissipation 
is also unimportant at this temperature. When Reu ^ 1> ohmic dissipation reduces the linear 
growth rate of the MRI (Sano & Miyama 1999) and the amplitude of the Maxwell stress at 
saturation (Sano & Inutsuka 2001). Note that we use a different definition of Reu than Fleming 
et al. (2000), which changes this critical value (see §5.2). The ratios shown in Figure [l] depend 



on the field strength, or c s /va (see eqs. |T3|] - |15|). The Hall parameter increases when c s /va 
increases because H/I = \X\/2 oc c s /va- For the case with c s /va = 1000, the Hall parameter is 
order unity, with a corresponding shift in the density range of the Hall regime to 10 10 <^n ri <^ 10 15 



cm 3 . 



For T = 1500 K, the Hall parameter \X\ ~ 1, and the Hall term can play an important role 
in the evolution of the MRI. The magnetic Reynolds number Reu = (O/I)^ 1 ~ 1 - 100 at this 
temperature. Therefore, both the Hall effect and ohmic dissipation are essential to the evolution 
of the MRI in dwarf nova disks when the temperature is less than 1500 K. 



2.2. Protoplanetary Disks 

The temperatures of protoplanetary disks are very low, so that the only sources of ionization 
are non-thermal, e.g., cosmic rays (Umebayashi & Nakano 1988; Gammie 1996), X-rays (Igea 
and Glassgold 1999), and the decay of radioactive elements. We shall briefly investigate which 
non-ideal MHD effects are important at different radii in the midplane of several disk models. 
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We assume a power-low distribution for the column density S and the temperature T, 

s « = s °(riu)"' (17> 

and 

T «= r »(riu)" 2 ' < 18 > 

where the mass of the central star is set to be M = M & . The minimum-mass solar nebula is 
chosen as a fiducial model; S = 1.7 x 10 3 g cm" 2 , p x = 3/2, T = 280 K, and p 2 = 1/2. The 
sound speed is c s (r) = [kT(r)/ fimn] 1 / 2 , where k is the Boltzmann constant and fi = 2.34 is the 
mean molecular weight. The angular velocity is tt(r) = (GMq/t -3 ) 1 / 2 , where G is the gravitational 
constant. For simplicity, cosmic rays and radioactive elements are considered as the only sources 
of ionization. The ionization rate due to cosmic rays, Ccr> is dominant because Cci? is about 
five orders of magnitude larger than the ionization rate due to radioactive elements Cr- But 
cosmic rays cannot penetrate into the disk more than the attenuation length xcr ~ 96 g cm -2 
(Umebayashi & Nakano 1981). Thus, the ionization rate at the midplane decreases dramatically if 
the column density of the disk exceeds XCR- The ionization rate at the midplane is given by 



CM ~ CcRexp 

where Ccr ~ 10~ 17 s^ 1 and Cr ~ 6.9 x 10~ 23 s" 1 



S(r) 

2XCR 



+ Cr , (19) 



The ionization fraction at the ionization-recombination equilibrium is approximately given by 



- = \ hr- ■ < 20 > 

n n V P n n 

where (3 = 1.1 x 10 _7 (T/300 K) _1 is the dissociative recombination rate (Millar, Farquhar, & 
Willacy 1997). The number density of electrons n e is a function of n n at the midplane, which is 

"»« = v^Sw) • (21> 

where H(r) = c s (r)/Q(r) is the scale height of the disk. For the magnetic field strength we simply 
assume c s /va = 1 or 10 at the midplane. 

Figure || shows the radial distribution of the ratios H/I, O/I, and A/ 1 at the midplane for 
the fiducial model with (a) c s /v\ = 1 and (b) c s /-i>a = 10. When c s /va = 1, the dominant term 
is the Hall effect within 10 AU, while ambipolar diffusion dominates in the outer parts of the 
disk. When the collision frequency jpi/Q = (A/I)~ l is greater than unity, the disk is linearly 
unstable for axisymmetric perturbations (Blaes & Balbus 1994). Thus the outer region of the disk 
is marginally unstable to the MRI. 

From equation ( |i~4| ) the ratio H/O is always larger than unity when the density of neutrals is 
n n <^ 10 18 (c s /ua) -2 cm -3 . Note that this is independent of the ionization fraction. Thus the Hall 
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effect is more important than ohmic dissipation almost everywhere in protoplanetary disks. The 
ratio H/O becomes larger as the density decreases, or the distance from the star r increases. The 
outer part of the disk has a higher ionization fraction because of its lower density. This means 
that both H and O decrease as r increases. We define ro and rjj as the radii where the ratios O/I 
and H/I equal unity. The Hall term is important within r < m- 

Protoplanetary disks can therefore be divided into three regions: (1) r > rjj, (2) ro < r < rjj, 
and (3) r < ro- In the outer part of the disk (r > rjj), both the Hall effect and ohmic dissipation 
are unimportant (\X\ < 1 and ReM 3> 1)- The critical radius rn is about 10 AU in the fiducial 
model with c s /vx = 1. At intermediate radii (ro < r < ru), the Hall term becomes important 
(\X\ > 1) but ohmic dissipation still can be ignored (ReM > 1)- For the region within ro, the 
ionization fraction is so low that both effects are important (\X\ 3> 1 and R&m ^ !)• I n the 
fiducial model the critical radius ro ~ 1 AU. In the innermost part of the disk r <^ 0.1 AU, thermal 
ionization is efficient and the gas becomes well coupled to the magnetic field. 

The critical radius rn is very sensitive to the field strength. For example, almost all of the 
disk within 100 AU is in the Hall regime for the case c s /va = 10, as shown by Figure ^b. The Hall 
parameter increases from 1 to 10 5 as the radius r decreases from 100 to 1 AU. The change in the 
critical radius ro is small. Within r = ro the magnetic Reynolds number is very small (1 - 10 ). 
The ratio A/ 1 is independent of c s /va- 

The dependence of r# and ro on the distribution of the column density is summarized in 
Table EL Disk models with pi = 3/2 and 1 are examined for two different disk masses Mdi s k- 
The dependence of the critical radii on p\ is very weak in these cases. As the disk mass is 
increased, the density at the midplane increases, so that the ionization fraction decreases and the 
critical radii shift outward in both cases. The density at the midplane also increases when the 
temperature decreases, because the scale height of the disk H is proportional to T 1 / 2 (see eq. 
plfl ). Furthermore the recombination rate increases as the temperature decreases. Thus a lower 
temperature means a lower ionization fraction. If a constant temperature T = 10 K is assumed 
everywhere in the disk, the critical radii become ro = 2.9 AU and rjj = 36 AU for the fiducial E 
distribution with c s /v\ = 1. 

In all of the models examined here, a large part of the disk (r < rjj ~ 10 — 100 AU) is in 
the Hall regime, while the inner regions (r < ro ~ 1 - 5 AU) are very resistive (ReM < 1). As 
discussed at the start of this section, this picture could be modified by the effect of dust grains 
(Wardle k Ng 1999; Sano et al. 2000). 



3. NUMERICAL METHOD 
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3.1. Equations and Algorithms 

We solve the equations of non-ideal MHD in a local Cartesian frame of reference corotating 
with the disk at the angular frequency O corresponding to a fiducial radius r$. The coordinates in 
this frame, written in terms of cylindrical coordinates (r, cp, z), are x = r — vq, y = ro((p — Fit), and 
z. The vertical component of gravity is ignored in this paper. For a small region Ar surrounding 
the fiducial radius tq, with Ar <C vq, the basic equations are given by 

^ + v • Vp = -pV • v , (22) 

dv VP J x B _ - 

+ v .\/ v = + 2f2xv + 2qti?xx , 23) 

at p cp 

de „ PV-v 4vrr,J 2 

+ „ . Ve = + —L— , (24) 

at p c z p 

— = V x [v x B : = V x £ , (25) 

at \ c en e ) 

where e is the specific internal energy and £ is the electromotive force (EMF). The term 2qQ 2 x in 
the equation of motion (]23| ) is the tidal expansion of the effective potential with a constant q = 3/2 
for a Keplerian disk. The gas is assumed to be the ideal, with pressure P = (7 — l)pe and 7 = 5/3. 
The induction equation (p5|) includes terms for both ohmic dissipation and the Hall effect. 

These equations are solved using a finite-difference code (Sano, Inutsuka, & Miyama 1999) 
in the local shearing box developed by Hawley, Gammie, & Balbus (1995). The hydrodynamics 
module of our scheme is based on a second order Godunov scheme (van Leer 1979) using a 
nonlinear Riemann solver that is modified to account for the effect of tangential magnetic fields. 
The evolution of the magnetic field is calculated with the Constrained Transport (CT) method 
(Evans & Hawley 1988), which guarantees the divergence free condition of the field, V • B = 0, 
is maintained. To extend the method to include non-ideal MHD effects, we use an operator 
split solution procedure. The MoC-CT technique (Stone & Norman 1992) is used to update the 
inductive term v x B in the EMF. The algorithm to update the Hall EMF, Sh = — J x B/en e , is 
described in detail in the Appendix. 



3.2. Initial Conditions 

Our simulations begin with Keplerian shear flow, so that in the corotating frame the azimuthal 
velocity is given by v y o = —qVtx. Two initial field configuration are considered in this paper: a 
uniform vertical field B z = Bq, and a zero-net flux vertical field B z (x) = Bq svl\{2'kx/ L x ), where 
L x is the box size in the radial direction. Except for the shear velocity, the initial state is uniform: 
p = Po and P = Pq . 
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The calculations are performed in a two-dimensional region in the radial-vertical (x-z) plane, 
in a volume bounded by x = ±H/2 and z = ±H/2, where H = (2/^) l l 2 c s /Q> is the scale height of 
the disk. Most of the runs use a standard grid resolution of 128 x 128 with uniform zoning. In the 
vertical direction, periodic boundary condition is used. For the radial boundary, a sheared periodic 
boundary condition (Hawley & Balbus 1992) is adopted. We choose normalizations with po = 1, 
H = 1, and = 10~ 3 . Then the sound velocity and gas pressure are initially 2c 2 /7 = 10 -6 and 
Pq = 5 x 10~ 7 , respectively. Initial perturbations are introduced as spatially uncorrelated pressure 
and velocity fluctuations. These fluctuations have a zero mean value with a maximum amplitude 
of \5P\/P = 10~ 2 and \Sv\/c s = 1(T 2 . 

This system is characterized by three non-dimensional parameters, /3q, Re mo, and Xq- The 
initial plasma beta (3q = {2 / ^)c 2 / v 2 A0 measures the initial field strength, where vao = Bq/(4ttp ) 1 / 2 . 
The magnetic Reynolds number ReMO = Vao/ 7 !^ se ^ s the value of the magnetic diffusivity r], which 
is assumed to be spatially and temporally constant in the calculation. The initial Hall parameter, 



*o = _: U " 2 , (26) 



cB Q 
2iren e0 v A0 

indicates the importance of the Hall effect and determines the initial number density of electrons 
n e Q. We assume the electron abundance is constant throughout the calculation and thus 
n e = n e Qp/pQ. Our simulations show that the density fluctuations in the nonlinear regime are not 
large, so that even if n e is assumed to be constant (i.e., n e = n e Q throughout the simulation) there 
is little difference to the results presented here. 



3.3. Linear Growth Rate of the Magnetorotational Instability 

To test the accuracy of our numerical algorithms, as well as to summarize the effect of the 
Hall term on the linear dispersion relation of the MRI, we first compare numerically measured 
growth rates with the predictions of linear theory. Balbus &; Terquem (2001) derive the dispersion 
relation for axisymmetric perturbations 5 tx exp(ikz + at) on a uniform vertical field B. Using the 
two non-dimensional parameters Reu and X defined by equations @ and (|l0|), the dispersion 
equation can be written as 



Reu ReM v ' 

k A k 2 



Co 



C 2 = R 2 + 2k 2 + -^ + '^-X[k 2 -A + k 2 X) , (28) 

k 2 -4 + k 2 (l + X)}+^- , (29) 

Re M 




where a = a/Q, k = kv\/£l, and k = k/Q, is the epicyclic frequency. A positive real root a implies 
an unstable exponential growth of the mode. 
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Figure ||a shows the dispersion relations for ReM = 100 and X = 4, 0, —2, and —4. Since 
R&M 3> 1) the effect on growth rates of the MRI is dominated by the Hall term. The sign of the 
Hall parameter X is determined by the direction of the magnetic field relative to ft. Positive 
(negative) X denotes that the field direction is in the same (opposite) sense of fl. Modes with a 
smaller wavenumber than the critical value k cr \t are unstable for the MRI. When X is positive 
the critical wavenumber decreases as the Hall parameter increases. The maximum growth rate is 
independent of X and has the value a ma x ~ 0.75, the same as the ideal MHD case. For X < 0, 
the critical wavenumber £; cr it — ► oo; the MRI has no characteristic scale in this case. On the other 
hand, for X < —4, no mode is linearly unstable to the axisymmetric MRI. 

The dispersion relations for ReM = 1 and X = 4, 0, —2, and —4 are shown in Figure ^o. 
Now the maximum growth rates are significantly reduced by ohmic dissipation. Fields oriented 
such that X is positive produce a larger growth rate than those with negative X. The critical 
wavenumber for negative X remains finite in this case, and X < —4 is always linearly stable for 
all axisymmetric modes. Thus the larger Hall parameter |X| enhances the growth rate of the MRI 
in some case (X ^ 1 at least at small k) and completely suppresses any growth in other case 
(X < —4), depending on the field geometry. 

Decreasing ReM further leads to even smaller growth rates and critical wavenumbers. For 
example, at ReM = 0.1 the maximum growth rate is <7 max = 0.11 (0.037) for the case of X = 2 
(—2), while the critical wavenumber is k cr n = 0.21 (0.12). These values are about an order of 
magnitude smaller than the ideal MHD case. 

We have measured the growth rate of the MRI using numerical simulations with uniform B z , 
Po = 400, and 256 x 256 grid zones. The simulations begin with very small random perturbations, 
|<5P|/Pinit = 10 -6 and |<5t>|/c s = 10 -6 , so that any growing modes should be well described by 
the linear analysis during the first a few orbits of evolution. A 2D Fourier decomposition of the 
radial velocity v x is carried out at frequent time intervals, and the growth rates are measured for 
modes with zero radial wavenumber, k x = 0, between orbits 2 and 2.5. These numerical growth 
rates are shown in Figure || by open circles. The numerical growth rates reproduce the analytical 
dispersion relation extremely well, each point is within a few percent of the predicted linear value. 
This provides a powerful validation of the numerical method. 

4. NONLINEAR EVOLUTION OF THE MAGNETOROTATIONAL 

INSTABILITY 

4.1. 2D Simulations with Uniform Vertical Fields 

In the ideal MHD case, the 2D MRI on a uniform B z evolves into an exponentially growing 
channel flow whose amplitude is unbounded (Hawley & Balbus 1992), although saturation at finite 
amplitude can occur when ohmic dissipation is included (Sano et al. 1998). We now examine how 
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these results are changed by inclusion of the Hall effect. 

The initial field geometry for all models in this subsection is a uniform vertical field B z = Bq, 
where -Bo is constant. Various models with different ReMO and Xq have been calculated and are 
listed in Table |2[ The critical wavelength A cr ;t = 2ir/k cr it and the maximum growth rate <r max are 
obtained from the dispersion equation (p7|). Hereafter (/) denotes the volume- averaged value of 
/, while ((/)) denotes the time- and volume-averaged value. The time averaging is taken over the 
last 20 orbits of the calculation. 

When the Hall parameter Xq is positive, the nonlinear evolution of the MRI is quite similar 
to the case Xq = 0. For all the models with ReMO <; 1 and Xq > 1, the nonlinear evolution of the 
MRI goes into the two-channel flow as in the ideal MHD cases. Figure |I| shows the time evolution 
of the magnetic field lines in the ReMO = 10 and Xq = 2 run (model Z5). The gray color contours 
denote the radial velocity v x normalized by the initial sound speed c s q. In this figure, time is 
given in orbits i ro t = 2n/fl. In the linear phase (until about 2 orbits), the modes with the most 
unstable wavelength Amri ~ 2ttva/Q dominate. However, because many other unstable modes are 
also growing, in the nonlinear regime they strongly interact leading to a transient phase of MHD 
turbulence. 

In the turbulent phase, the vertical component of the field is amplified, and the characteristic 
length of the MRI in the vertical direction becomes longer. Thus, at late times the magnetic 
field and velocity are dominated by larger structures. The final state shown by the last panel of 
Figure || is quite similar to the ideal MHD case (Hawley & Balbus 1992). The magnetic field 
lines are almost horizontal, and two radial streams, one inward and the other outward, occupy 
the computational box. The increase of the magnetic energy continues without saturation until it 
exceeds 10Po and the calculation is stopped. 

Figure || shows the time evolution of the volume-averaged magnetic energy (B 2 /87r) for 
models with Xq = 2, 0, and —2 (model Z5, Z6, and Z7). The magnetic Reynolds number of these 
models is ReMO = 10. As mentioned above, the Xq = 2 and runs show unbounded exponential 
growth. Interestingly, the Xq = —2 run shows amplification of the magnetic energy by an order of 
magnitude during the linear phase, followed by saturation at finite amplitude. This amplitude is 
sustained for at least 50 orbits. 

We find that the growth of the MRI saturates for all models with Xq < 0. Figure || shows 
images of typical time evolution of the field lines in a saturated model (Z7). The angular 
momentum perturbations 5v y = v y — v v q are also shown by gray color contours. Obviously, smaller 
scale fluctuations are amplified in the linear phase compared to the model with Xq = 2 (shown 
in Fig. |j). Larger scale structures do not dominate at the nonlinear stage for Xq < 0, because 
small wavelength perturbations are always unstable even after the field is amplified by the MRI 
(see §5.1). 

The nonlinear stage of model Z7 is MHD turbulence as shown by the last panel of Figure ||. 
The efficiency of angular momentum transport is measured by the turbulent stress w xy , which is 
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the sum of the Maxwell stress, —B x B y /4TT, and the Reynolds stress, pv x 5v y . The a parameter 
of Shakura & Sunyaev (1973) can then be expressed as a = w xy /Po. The Maxwell stress is 
well-correlated with the magnetic energy, so that the time evlution of the stress is similar to that 
of the magnetic energy shown by Figure We find that the saturation level of the Maxwell and 
Reynolds stresses are almost the same in model Z7. The time- and volume-averaged stresses 
are ((-B x By/4ir))/P = 1.4 x 1(T 3 and ((pv x 5v y )) / P = 1.3 x 1(T 3 . The turbulence in this 
case is slightly anisotropic: the components of the magnetic and perturbed kinetic energy are 
{{Bl)) = 4.6((B 2 X )) = l.8((B 2 z )) and ((pSv 2 y )) = 1.2«» = U«. 

The effective Hall parameter in the nonlinear regime is defined as 

XcS = cBzQ = 2cpBzn (30) 
c 2Tren e v\ en e B 2 

Because the Hall parameter is inversely proportional to the field strength, the volume-averaged 
value |(X e ff)| decreases as the field is amplified. The time-averaged value is ({X e g}} = —0.16 for 
model Z7, which started with Xq = —2. At the end of model Z5, the volume-averaged effective 
Hall parameter is (X e g) ~ 10 -4 . The sign of (X c s) remains unchanged for both models. Thus 
the importance of the Hall term becomes smaller as the MRI amplifies the field strength in the 
disk. When the field is amplified, the magnetic Reynolds number at the nonlinear stage (v\ z /rjQ) 
increases, and then ohmic dissipation also becomes less effective (Sano et al. 1998). 

When ReMo < 1, the growth of the MRI saturates even without the Hall term (Sano et 
al. 1998). Thus, we find the Hall term has less effect in this case: the nonlinear evolution with 
Xq = ±2 is quite similar to that of Xq = 0. Ohmic dissipation controls the MRI in these cases, 
and the dependence on Xq is very small except for the linear growth rate. Figure [7| shows the 
time evolution of the azimuthal component of the magnetic energy (B y /8ir) for the models with 
R&M0 = 0.1 and Xq = 2, 0, and —2 (model Z13, Z14, and Z15). The linear growth phase continues 
until a few tens of orbits because of the small growth rates in these models. The magnetic energy 
is amplified up to the initial thermal energy ~ Pq for all the models. However the energy is not 
sustained, and finally dies out after 100 orbits. The Maxwell stress at the end of the calculation is 
very small for these cases, ((— B x B y / 4tt)) / Pq ~ 10~ 9 . 

Although the magnetic field returns to its initial configuration after 100 orbits, we do not see 
the MRI re-emerge. In 3D, re-emergence of the MRI on the background field (which can never 
be dissipated) leads to strong fluctuations in, e.g., the Maxwell stress (Fleming et al. 2000). In 
these 2D simulations, the magnetic and perturbed kinetic energy is nearly comparable at the end 
of the linear growth phase, so that the kinetic energy is {p8v 2 /2) ~ 0.2P > £o/8tt. After that, 
the magnetic energy starts decaying while the kinetic energy remains larger than its initial value. 
The kinetic energy is contained mostly in velocity along the field (vertical) direction. This flow is 
subthermal but super Alfvenic; this may suppress further growth of the MRI. 

The evolution of the MRI on uniform vertical fields in 2D is summarized by Figure B The 
character of the nonlinear regime in the 2D parameter space defined by ReMo and Xq is denoted 
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by open circles for simulations that result in unbounded growth of the channel solution, filled 
circles for simulations that result in saturation and eventual decay, and crosses for simulations 
that show no growth at all. We find that when ReMo <; 1 and Xq > 0, an unbounded two-channel 
flow emerges in the nonlinear regime. Saturation at finite amplitude occurs when Xq < or 
R&M0 ^ 1- Finally, if Xq < —4, all modes are stable: no linear growth can be seen for the models 
with Xq = —5 (model Z4, Z8, Z12, and Z16) shown by crosses in the figure. 

4.2. 2D Simulations with Zero-Net Flux Vertical Fields 

Next we consider models with a zero-net flux vertical field, B z (x) = Bq sin(2-7rx/L x ). In the 
ideal MHD case with this initial field, the MRI evolves into a transient phase of MHD turbulence 
that eventually dies away (Hawley & Balbus 1992). For this initial field, the Hall parameter is 
given by X(x) = Xq/ sm(2nx/ L x ) where Xq = 2cpQQ/en e QBQ. Thus, the region x < has positive 
X while the region x > has negative X, and the minimum of is Xq. Table || lists the 

models computed with a zero- net flux. The initial plasma beta (3q, the critical wavelength A cr i t , 
and the maximum growth rate <7 max are given for B z = Bq and X = Xq. 

Figure ^ shows images of the magnetic energy in model S2, with ReMo = 100 and Xq = 2. 
The velocity field is also shown by arrows. Most of the region with negative X is stable because 
the Hall parameter X < —4. The growth of the MRI is seen mainly in the region with positive X, 
as evident in the first panel of Figure ||. The growth in x < gradually affects the structure of the 
region x > 0. The amplified magnetic field spreads across the entire computational domain by 3 
orbits, after that the amplified magnetic energy is sustained for at least 100 orbits. 

Figure [H] shows the time evolution of the magnetic energy for ReMo = 100, 10, and 1 with 
and without the Hall term (Xq = and 2). For the ReMo = 100 and 10 runs, the magnetic 
energy is sustained until 100 orbits. In the turbulent state the magnetic and the perturbed kinetic 
energy are equi-partitioned. Since the vertical component of the magnetic energy, which is the 
seed of the MRI, continues decaying throughout the evolution, the activity of the turbulence 
gradually weakens over time in both runs. Therefore, active MHD turbulence is a transient 
phenomenon in these cases. When ReMo = 100, the Maxwell stress at the end of the calculation 
is ((-B x B y /4Tr))/Po = 1.6 x 10~ 5 and 2.4 x 10~ 8 in the X = and 2 run, respectively. Only the 
toroidal field remains at late phases of the evolution. In 2D, a disk with a purely toroidal field is 
stable to the MRI. However, non-axisymmetric unstable modes could grow from the toroidal field 
and sustain MHD turbulence in 3D. 

The evolution of more resistive models with Re mo = 1 and Xq = and 2 is quite different. 
The magnetic energy is amplified for several orbits, but then dies out. The decay time of the 
magnetic energy is about 5 orbits. The diffusion time tdiff of magnetic fields for a length scale I is 
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given by 

Thus the decay time corresponds to the diffusion time tdiff with I ~ 0.1H ~ A cr it . 



The Hall parameter can take values up to about 100/i?eA/o (see eq. [fi6| ) when the disk is 
in the Hall regime. Thus, a small magnetic Reynolds number ReMo may give a very large Hall 
parameter Xq. Since the growth rate of the MRI is enhanced by large Xq at small Re mo, it is 
important to examine whether significant growth can occur in models with very large Xq. Models 
S7 and S8 have Rcmo = 1, and Xq = 10 and 100 respectively. The evolution of the Xq = 10 run 
shows decay of the magnetic energy similar to the Xq = 2 run (model S6). On the other hand, 
for Xq = 100, active turbulence is initiated and sustained for at least 10 orbits. Figure 11 shows 
the time evolution of each component of the magnetic field in this run. As in the other active 
turbulent models, the Maxwell stress is strongly correlated with the magnetic energy and the time 
evolution of the stress is resemble to that of the magnetic energy. The stress is amplified up to 
10~ 3 i-b an d sustained until 10 orbits. This saturation level is comparable to or slightly larger than 
those in the less resistive models (Re^o = 100, 10) at the same phase. 

For more resistive models with Re mo = 0.1, the nonlinear evolution in all cases (models S9 
- S12) shows rapid decay of the magnetic field even with very large Xq. The decay timescale of 
the magnetic energy is very short, about an orbit. In the Xq = 1000 run with Rcmo = 0.1 (model 
S12), the critical wavelength of the MRI is longer than the scale height of the disk H, or the 
vertical box size. No linear growth of the MRI can be seen in the calculation, nor would the MRI 
be present for these parameters in actual disks. 



5. DISCUSSION 

5.1. Interpretation of the 2D Evolution 

The axisymmetric MRI evolves into a channel flow or MHD turbulence, depending on the 
parameters Re mo and Xq (provided linearly unstable modes exist at all). The results summarized 
in Figure || reflect the linear properties of the MRI. From the dispersion relation, the parameter 
space Xq > —4 can be divided into three regions, as shown by dotted lines in Figure ||. When 
R&M0 ^ 1 an d Xq > 0, the characteristic wavelength of the MRI is proportional to va/Q, and 
a critical wavelength exists. For this case, the evolution of the MRI shows an inverse cascade of 
the magnetic energy, and a two-channel flow emerges without saturation. This is because the 
characteristic scale increases as the field strength is amplified by the growth of the MRI. In the 
second region (ReMo ^ 1); the dispersion relation has a critical wavelength, but the most unstable 
wavelength is proportional to tj/va, that is inversely proportional to the field strength. In the 
third region (Xq < 0), there is no characteristic scale of the MRI. Then smaller scale fluctuations 
are always unstable. Therefore, models in these last two regions show no evidence for an inverse 



-16- 



cascade; instead they evolve into MHD turbulence. 

Models with (5 = 800 also show the same characteristics as the models plotted in Figure || 
implying these results are independent of the initial field strength. Estimation of the saturated 
amplitude of the Maxwell stress, which determines the efficiency of angular momentum transport 
in actual accretion disks, requires 3D simulations of the MRI. In 3D simulations that include only 
the ohmic dissipation, the saturation amplitude of the Maxwell stress is larger if the MRI has 
an inverse cascade (Sano & Inutsuka 2001). Local 3D simulations of the MRI in Hall MHD are 
presented in Sano & Stone (2002). 

5.2. Definition of Magnetic Reynolds Number 

In this paper we have defined the magnetic Reynolds number as Reu = v a/v^i an d used 
ReMo to denote the value of this number in the initial state. Both a linear analysis of the MRI 
including ohmic dissipation, as well as local 3D simulations with an initial uniform vertical field 
reveal that the critical number required to generate significant MHD turbulence and angular 
momentum transport is ReMo ~ 1 (Sano & Inutsuka 2001). If ReMo > 1 the evolution and 
saturated state of the MRI is little changed from the ideal MHD case, whereas if ReMo < 1 the 
growth rates and amplitude of the saturated state are both significantly reduced compared to the 
ideal MHD case. 

In actual accretion disks, there is a minimum value for Reu that results from the requirement 
that the critical wavelength of unstable modes in the very resistive regime Amri ~ v/ v A be less 
than the disk thickness H (Sano & Miyama 1999). This requires Re mo > va/c s , which can also 
be used as a stability criterion (Gammie 1996; Igea & Glassgold 1999; Sano et al. 2000). If the 
field strength of accretion disks is subthermal and va/c s < ReMo < 1) the MRI will be present, 
but the growth rate, saturated amplitude, and angular momentum transport rate will all be 
reduced compared to the ideal MHD case. In practice, the largest value for the ratio va/c s that 
can be reached in an accretion disk is probably unity (otherwise the field will escape the disk 
via buoyancy). In this case, the magnetic Reynolds number becomes Re' M = c^/r]Q (Gammie &; 
Menou 1998; Fleming et al. 2000), with the critical value measured from 3D simulations about 
10 4 . This last form is independent of the magnetic field strength in the disk, and therefore can be 
measured more easily with observations. However, even if the observed Re' M is below the critical 
value 10 4 , the MRI may still operate if 0.01 <^va/c s < 1. 

6. SUMMARY 

We have shown that in both dwarf nova disks in quiescence, and in protoplanetary disks 
around young stars, the Hall term can be important to the linear properties of the MRI, and 
therefore must be included in realistic models of these systems. Next, we investigated the nonlinear 
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evolution of the MRI in Hall MHD using axisymmetric 2D simulations. Although 3D simulations 
are essential for evaluating the efficiency of angular momentum transport in accretion disks, 
understanding the 2D evolution will be very helpful to the interpretation of 3D results. Fully 3D 
simulations of the nonlinear regime of the MRI in Hall MHD are presented in a companion paper 
(Sano & Stone 2002). 

Our findings are summarized as follows. 

1. For dwarf nova disks in quiescence, both the Hall effect and ohmic dissipation are very 
important when the temperature is less than about 1500 K. The importance of these effects 
is very sensitive to the temperature of the disk. 

2. In the region r <^ m in protoplanetary disks, the Hall term dominates over ohmic dissipation 
and ambipolar diffusion, assuming dust grains have settled out to the midplane. The critical 
radius r# ~ 10 - 100 AU, and depends mainly on the magnetic field strength in the disk. 
Inside a radius of ro ~ 1 - 5 AU, ohmic dissipation suppresses all modes of the MRI, except 
within about 0.1 AU where thermal ionization becomes important. If dust grains with the 
same size distribution as interstellar grains are well mixed throughout the disk (i.e., they 
have not settled to the midplane) , ro increases dramatically, and the entire disk is dominated 
by ohmic dissipation. 

3. For models with an initially uniform B z , the MRI evolves into either a two-channel flow 
which grows without bound (for Reuo ^ 1 and Xq > 0), or saturates as MHD turbulence 
(for ReMo 1 or Xq < 0). If Reuo 1, saturation of the MRI occurs, but the amplified 
magnetic field eventually dies out. 

4. For models with an initial zero-net flux vertical field, the MRI saturates as MHD turbulence, 
which either is sustained or eventually dies away depending on ReMo- In this case, the 
evolution is determined mostly by the effect of ohmic dissipation and has little dependence 
on the Hall parameter Xq. 

We thank Steven Balbus, Caroline Terquem, and Neal Turner for helpful discussions. 
Computations were carried out on VPP5000 at the National Astronomical Observatory of Japan 
and VPP700 at the Subaru Telescope, NAOJ. 

A. Algorithm for the Calculation of the Hall Term 

We adopt an operator split solution procedure for the update of the magnetic field. Following 
the constraint transport method (Evans & Hawley 1988), the Hall electromotive force (EMF) 
Eh is calculated at the zone boundaries. For simplicity the ID algorithm (in the x direction) is 
described here; the extension to multi-dimensions is straightforward. In the ID case, the tangential 
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components of the magnetic field and EMF are denned at the zone center and zone boundary, 
respectively, and the x component of the magnetic field is constant B x . 

The change of the magnetic field due to the Hall EMF is given by 

— = X7 xS H = V x [-Q H (VxB)xB] , (Al) 



where 



where Qh = c/47ren e is assume to be constant. Each component of the induction equation (Al) is 
written as 

dBy_ ae BM _ &b„ 

dt dx x dx 2 ' 

dB z _ d£ Hy _ 8 2 B y 

~df - "aT - - QhBx -^ ' (A3) 

BE 

£ Hy = -QhB x -^ , (A4) 
BE 

£hz = -QhB x ^ . (A5) 

For the update from t n to t n+1 , the partially time-advanced values of the Hall EMF S 1 ^ 1 ^ 2 at 
t n+l l 2 are needed. From equations (A4) and (|A5|), the time-advanced Hall EMF at a zone x = i 
and t n+l / 2 is given by 

Tjii+l/2 _ R n+l/2 

£H + yf = Qulh A; , ' , (A6) 

R n+l/2 _ R n+l/2 

b Hz,i ~ ~QhB x — , (A7) 

where Ax is the grid scale. To evaluate B n+l l 2 , we use equations (A2) and (A3); the induction 
equation including only the Hall EMF. Then the partial update of B y and B z is given by 

R n+l/2 _ Rn At B z,i+1 ~ 2B z,i + B z,i-1 ,.„x 
B y,i ~ B y,i + -^QhSx ^2 ' ( A8 ) 



—Qh b x ^ 

where At is the time step. 



R n+l/2 _ R n /OR y ' i+1 M y ^~ l ( \d\ 

B z,i - B z,i ~ ~itQhB x — — , (A9) 



Since our update of the Hall term is time explicit, we also must add a new time step 
constraint. The frequency of the whistler wave is approximately given by u ~ QnB x k 2 , where k is 
the wavenumber, thus the group velocity is 2QuB x k. The new constraint on the time step is then 



Ax (Ax) 2 
At ~ 2Q^B~k = 4^B~ X ' (A10) 



where we assume k = 2tt/ Ax. 
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In addition to checking the numerical growth rates of the MRI in Hall MHD with the 
predictions of linear theory (Fig. |3|), we have also tested our algorithm by comparing the 
numerically measured phase velocity of whistler waves with the linear dispersion relation. For 
this test, we adopt a uniform density p and a uniform magnetic field parallel to the x axis, B x 
{By = B z = 0), and consider the propagation of linear waves proportional to expi(kx — ut) in 
the incompressible limit. The dispersion equation without the Hall effect is uj 2 = v\k 2 where 
v\ = B 2 /4irp. If the Hall term is included, the dispersion equation is then written as 

(u 2 -v 2 A k 2 ) 2 = Q 2 H B 2 kW . (All) 



Figure |12| shows the dispersion relation obtained from equation ( All ). Due to the Hall effect, the 
Alfven wave {u/vxk = 1) splits into left- and right-circularly polarized waves. The right (whistler) 
wave can propagate much faster than the Alfven wave. Numerically obtained phase velocities are 
also shown in the figure by open circles. Each mode is calculated using a ID simulation of the 
propagation of a linear wave with 32 grid zones per a wavelength. The numerically measured 
phased velocities agree with the dispersion relation very well (better than 1%). This shows that 
our algorithm accurately captures the characteristics of Hall MHD. 
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Table 1. Critical Radii in Models of Protoplanetary Disk 1 



Cs/va = 1 c s /v A = 10 



So 




Pi 


M disk [M ] b r c 


, [AU] 


r H 


[AU] 


r [AU] r H 


[AU] 


1.7 x 10 3 




3/2 




0.024 




1.4 




8.9 


4.0 


76 


1.7 x 10 4 




3/2 




0.24 




4.7 




16 


9.7 


80 


3.3 x 10 2 




1 




0.024 




0.58 




8.4 


2.9 


76 


3.3 x 10 3 




1 




0.24 




2.5 




14 


7.4 


84 


a The temperature distribution is 


assumed to be Tq 


= 280 K and p 2 = 


1/2 for all models. 




b The disk mass -M^k 


is integrated from 0.1 to 100 AU. 
















Table 2. Uniform B z 


Simulations 






Model 


Size 


Grid 


A) 


ReMo 


^0 


Acrit/ H 


o"max/^ Orbits 


Zl 


H 


x H 


128 


x 128 


3200 


100 


2 


0.11 


0.75 


4.5 


Z2 


H 


x H 


128 


x 128 


3200 


100 





0.064 


0.74 


5.9 


Z3 


H 


x H 


128 


x 128 


3200 


100 


-2 


0.0 


0.70 


50 


Z4 


H 


x H 


128 


x 128 


3200 


100 


-5 






50 


Z5 


H 


x H 


128 


x 128 


3200 


10 


2 


0.11 


0.72 


5.0 


Z6 


H 


x H 


128 


x 128 


3200 


10 





0.064 


0.70 


5.4 


Z7 


H 


x H 


128 


x 128 


3200 


10 


-2 


0.0 


0.60 


50 


Z8 


H 


x H 


128 


x 128 


3200 


10 


-5 






50 


Z9 


H 


x H 


128 


x 128 


3200 


1 


2 


0.12 


0.51 


6.4 


Z10 


H 


x H 


128 


x 128 


3200 


1 





0.091 


0.43 


9.7 


Zll 


H 


x F 


128 


x 128 


3200 


1 


-2 


0.064 


0.28 


50 


Z12 


H 


x F 


128 


x 128 


3200 


1 


-5 






50 


Z13 


2H 


x 2H 


128 


x 128 


3200 


0.1 


2 


0.54 


0.11 


100 


Z14 


2H 


x 2# 


128 


x 128 


3200 


0.1 





0.64 


0.074 


100 


Z15 


2H 


x 2H 


128 


x 128 


3200 


0.1 


-2 


0.91 


0.037 


100 


Z16 


2H 


x 2# 


128 


x 128 


3200 


0.1 


-5 






100 
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Table 3. Zero-Net Flux B z Simulations 



Model 


Size 


Grid 


00 


ReMo 


^0 




"'max/^ 


Orbits 


SI 


H x 


H 


128 


x 


128 


3200 


100 





0.064 


0.74 


100 


S2 


H x 


H 


128 


x 


128 


3200 


100 


2 


0.11 


0.75 


100 


S3 


H x 


H 


128 


X 


128 


3200 


10 





0.064 


0.70 


100 


S4 


H x 


H 


128 


X 


128 


3200 


10 


2 


0.11 


0.72 


100 


S5 


H x 


H 


128 


X 


128 


3200 


1 





0.091 


0.43 


100 


S6 


H x 


H 


128 


X 


128 


3200 


1 


2 


0.12 


0.51 


100 


S7 


H x 


H 


128 


X 


128 


3200 


1 


10 


0.22 
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Fig. 1. — Ratios of non-ideal MHD terms to the inductive term given by equations (j9|) through 
(11) for the case of a dwarf nova disk in quiescence for (a) T = 3000 K and (b) T = 1500 K. 
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Fig. 2. — Ratios of non-ideal MHD terms to the inductive term given by equations (P) through 



(11) for the case of a protoplanetary disk with (a) c s /va = 1 and (b) c s /v^ = 10. 
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Fig. 3. — Dispersion relation for axisymmetric MRI modes for the case with (a) Rem = 100 and 
ReM = 1- Open circles depict the numerical growth rates measured in 2D numerical simulations. 
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Fig. 4. — Magnetic field lines (solid) and radial velocity, 6v x /c s q (gray color) in model Z5 ((3q = 3200, 
Reuo = 10, and Xq = 2) at orbits 1.5, 2, 3, and 5. 




Fig. 5. — Time evolution of the volume-averaged magnetic energy (B 2 /Stt)/Pq for model Z5 
(Xq = 2), Z6 (Xq = 0), and Z7 (Xq = —2). The plasma beta and the magnetic Reynolds number 
of these models are (5q = 3200 and ReMo = 10. 
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Fig. 6. — Magnetic field lines (solid) and angular momentum perturbation, 5v y /c s o (gray color) in 
model Z5 ((3 = 3200, Re M o = 10, and X = -2) at orbits 2, 3, 5 and 10. 




Fig. 7. — Time evolution of the azimuthal component of the magnetic energy (By/87t)/Po for model 
Z13 (Xq = 2), Z14 (Xq = 0), and Z15 (Xq = —2). The plasma beta and the magnetic Reynolds 
number of these models are {3q = 3200 and ReMo = 0.1. 
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Fig. 8. — Final state for calculated disk models with a uniform vertical field. Filled circles denote 
models in which the growth of the MRI saturates as MHD turbulence in the nonlinear stage. Open 
circles denote models which evolve into a channel flow that grows without bound. Crosses denote 
models which are linearly stable to the MRI. The parameter space (Xq, ReMo) is divided into four 
regions shown by the dotted lines according to the characteristic of the linear dispersion relation 
(see §5.1). 
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Fig. 9. — Magnetic field energy, log(i? 2 /87rPo) (9 ra D color) and velocity field (arrows) in model S2 
(/3o = 3200, Re M o = 100, and X = 2) at orbits 1.5, 2, 3, and 5. 
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Fig. 10. — Time evolution of the magnetic energy ((B 2 /8tt))/Pq for model SI (ReMo = 100 and 
X Q = 0), S2 (Re M o = 100 and X = 2), S3 (Re mo = 10 and X = 0), S4 (i?e M o = 10 and X = 2), 
S5 (ReMo = 1 and Xo = 0), and S6 (Reuo = 1 an d = 2). The plasma beta of these models is 
A) = 3200. 
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Fig. 11. — Time evolution of each component of the magnetic energy and the Maxwell stress in 
model S8 (/% = 3200, Re M o = 1, and X = 100). 




k 

Fig. 12. — Dispersion relation of the left- and right-circularly polarized wave traveling along field 
lines. The Alfven wave for the case without the Hall effect (uo/v^k = 1) is shown by dotted line for 
comparison. Open circles denote numerically measured phase velocity in ID simulations. 
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